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We calculate the low-temperature thermodynamic quantities (magnetization, correlation func- 
tions, transverse and longitudinal correlation lengths, spin susceptibility, and specific heat) of the 
frustrated one-dimensional spin-half J1-J2 Heisenberg ferromagnet, i.e. for J2 < 0.25|Ji|, in an ex- 
ternal magnetic field using a second-order Green-function formalism and full diagonalization of finite 
systems. We determine power-law relations for the field dependence of the position and the height 
of the maximum of the uniform susceptibility. Considering the specific heat at at low magnetic 
fields, two maxima in its temperature dependence are found. 
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I. INTRODUCTION 

Frustrated low-dimensional spin systems have at- 
tracted increasing attention*^ They represent an ideal 
playground to study the influence of strong (thermal 
and quantum) fluctuations on thermodynamic quanti- 
ties. The one-dimensional (ID) J1 — J2 Heisenberg model 
with ferromagnetic nearest-neighbor (NN) coupling Ji 
and antifcrromagnetic next-nearest-neighbor (NNN) cou- 
pling J 2 has been studied extensively over the last years, 
see, e.g., Refs. HUTU It has been pointed out that this 
model is an appropriate starting point to describe ex- 
perimental results for the family of the quasi-lD edge- 
shared chain cuprates, such as LiVCuVi, L1CU2O2, 
NaCu 2 2 , Li 2 ZrCuOi, and Li 2 Cu0 2 ^^ The corre- 
sponding Hamiltonian of this system with an external 
magnetic field is given by 



H = J l Y J S i S j + J 2 SiSj - h $i 



(1) 



(id) 



[id] 



where runs over the NN and over the NNN 

bonds. For the exchange constants we assume J\ < 
and J 2 > 0. In most of these materials quite large val- 
ues of J2 are realized leading to incommensurate spiral 
in-chain spin correlations at low temperatures. Magnetic 
field effects can be quite strong j£ & 12 : 26 ~ — in particular, 
if J2 is near the quantum critical point Jf = |Ji|/4, 
where the transition between the ferromagnetic and in- 
commensurate spiral ground state takes place. On the 
other hand, some materials considered as ID s = 1/2 
ferromagnets, such as the copper salt TMCuC ) 29 : 30 the 
organic magnets p-NPNN 3 ^ and /3-BBDTA-GaBr 4 ^ 
may have a weak frustrating NNN exchange interaction. 
Moreover, in Refs. and [33 L12CUO2 and L%2ZrCuOi 
were identified as quasi-lD s = 1/2 frustrated ferromag- 
nets with J 2 ~ 0.2|Ji|. Although for J 2 < Jf the ground 
state of the model ([T]) is ferromagnetic, it has been shown 
recentl y 10 : 36 that the low-temperature thermodynamics 



is str ong ly influenced by the frustrating J2 ■ 

~ 1 it has been found that for unfrustrated ID 



Moreover, in 

Refs.|37|and|3£ 

quantum ferromagnets at weak magnetic fields the spe- 
cific heat exhibits an additional low-temperature maxi- 
mum indicating a separation of two energy scales. In- 
terestingly, this field-induced low-temperature maximum 
appears only in ID systems and for the small spin quan- 
tum numbers, s = 1/2 and s = L 37 : 38 Therefore, it can 
be considered characteristic feature of ID quantum 
ferromagnets. 

In this paper we consider the interplay between frus- 
tration and magnetic field in ID s = 1/2 ferromagnets, 
i.e. we study the model (Ji} in the parameter regime 
J 2 < |Ji|/4, where the ferromagnetic ground state is 
realized. In difference to previous investigations) 37 : 38 
neither the Bethe ansatz nor the quantum Monte Carlo 
method can be used. Hence, we employ (i) a second- 
order Green-function method (GFM) for infinite chains 
and (ii) exact diagonalization (ED) of finite chains with 
A = 16 and N = 20 spins imposing periodic boundary 
conditions to calculate thermodynamic properties. It has 
been shown in Refs. H0II37I and[38| that both techniques 
may provide reliable results for the problem under 
consideration. 

The Green-function technique for Heisenberg spin 
systems used here was initially introduced by Kondo 
and Yamaij in a rotation-invariant formulation^ It 
has been further developed and successfully applied 
to low-dimensional quantum systems over the last 
decadei 10 i 37 ' 38 i 40 : 41 The specific version of the method 
appropriate for spin systems in a magnetic field was de- 
veloped in Ref. [3|| 

The paper is organized as follows. In Sec. |TT] we il- 
lustrate the main features of the GFM. In Sec. Mil we 
present our results for the magnetization, the suscepti- 
bility, the spin-spin correlation functions, the correlation 
length, and the specific heat. A summary of our results 
is given in Sec. IIV1 
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II. SECOND-ORDER GREEN-FUNCTION 
THEORY 

To calculate the longitudinal and transverse spin cor- 
relation functions and other thermodynamic quantities, 
we use two-time retarded commutator Green functions^ 
First we determine the longitudinal spin correlation func- 
tions from the Green function ((S z ; SZ q ))ui = — Xq Z ( w )> 
where Xq Z is the longitudinal dynamic spin suscepti- 
bility. The equation of motion reads 



co 2 ((S z ;S z _ q ))^ = M zz + ((-S z ;S z _ q }}, 



(2) 



M zz 

q 



-2 J n C-+ (1 - cosnq) , (3) 



where C n + = (S S£). To approximate the second 
deriv ative — S z we use a decoupling scheme proposed in 



(4) 



Refs. I37H41I which reads 



S Z S+S~ =a zz k (S+S~)S z 



where for the vertex parameters a z %. we assume a z %, 
af z if j, k are NN, and a^f. 

-S z = «) 2 « and 



l\ z otherwise. We obtain 



M z 



,zz\2 



(5) 



Among others, Eq. ((9]) can be used to determine the ver- 
tex parameters. 

To calculate the transverse correlation functions, we 
use the Green function ((S q ; SZ q ))u = ~Xq~ where 
Xq~ (w) is the transverse dynamic spin susceptibility. 
Here the equations of motion read 



uj((S+;S_ q )) b 



2(S*) + ((iS+;SZ q )) u , (11) 
M+- + u((-S+;Sl q )) u . (12) 



The moment 



} is given by 



Ml 



-4 J n C n (l- cosnq) + 2h(S z ), (13) 



n=l,2 



where C n = \C n + + C zz . In —S+ we decouple the prod- 
ucts of three operators as 

S-S+S+ = afj (SrSf)S+ + a+ fc " (SZS+)S+, (14) 
where again for the vertex parameters af~ we assume 



"Jfc 

We obtain 



with 



j if j, k are NN, and a^Z~ = a\ otherwise 



s: 



[(u+-Y-h 2 ]S+ + 2hi& q 



(15) 



with 



«) 2 = 



^2 JnJmO- ~ cosnq) (6) 

n,m(=l,2) 

x [#"m + 4< Z C-+(1 - cosmg)] , 



where Jf£ n - 1 + 2(af C7+ - 3o^C~+), i^ z 2 
2(af C 3 " + - afCf +), and A'f^ = iff 2 + 4(afCf + 



a| z C^~ ). Applying the spectral theorem^ to the Green 
function |J5J), the longitudinal correlation functions of ar- 
bitrary range n, C zz = (SqS z ), can be calculated b y 37 ' 38 



q^Q 

M zz 

C fHSV^^[l + 2n«)] 



(7) 
(8) 



where n(ui zz ) = (e u i l T — lj is the Bose function. 

For n = 0, we get the sum rule C§ z = 1/4. As shown in 
Ref. [HI the isothermal and the uniform static Kubo sus- 
ceptibility Xq z agree at arbitrary fields and temperatures. 
Using Eqs. (J3j> , (0, and (|6|) we get the relation 



d(S z ) 
dh 



n=l,2 
n,m(=l,2) 



(9) 
(10) 



("t ) 2 = E J n Jm(l - cosnq) (16) 

n,m(=l,2) 

x Km + 4a+~C7 n (l - cosmg)] , 



where = 1 + 2(a+-C 2 „ - 3a+-C„), tf+J = 

2{a+-C 3 - af~Ci), and if+f = K+Z + 4(a+ J Ci - 
a\ C 2 ). Finally, the resulting Green functions are 



((S q ; S_ q )) b 
((iS£; S_ q )) lJ 



t=l,2 

E 



.4 



where 



J <?1,2 



A 9l , 2 = (S z ) ± -V- (M+- - 2h(S z )) 



(17) 
(18) 



(19) 
(20) 



As argued in Ref. \3aL a divergence in the transverse static 
spin susceptibility Xq ( w = 0), which signals a phase 
transition, could appear if w 9c2 = 0, i.e., ui q ~ = h. Since 
for nonzero fields, the Heisenberg ferromagnet (JXJ) docs 
not describe a phase transition, xj ( w = 0) has to be 
finite at all q. To ensure this, we require A Qc 2 = which 
results in the regularity condition 



h{S z ) = -2j2 J " C » (1 - cosnq c ) 



(21) 



71=1,2 



3 




0.1 0.2 0.3 0.4 0.5 
T 

FIG. 1: Magnetization for J2 = 0,0.1,0.15 (from top to bot- 
tom) and h = 0.005,0.03,0.1 versus temperature [GFM - 
lines; ED (N = 20) - open symbols]. 

Following Rcfs. |38| and |4l|, we assume this condition to 
be valid also for larger fields, where h > us^ for all q, to 
guarantee the continuity of all quantities. 
Applying the spectral theorem to the Green functions 
([T7|) and (fl8| . we obtain the transverse correlation func- 
tions C~+ = {SoS+) = iJ2 q C~ + e^ and C~+ = 

(SqIS^) = jj J2 q C~ + e lqn of arbitrary range n with the 
structure factors 

C~+ = A qi n(oj qi ),C- + = ^2 u q iA qi n(u}qi). (22) 

i=l,2 i=l,2 

From the s = 1/2 operator identity S~ Sf~ + S z = | we 
find the sum rule 

<<?*) = 1-C -+. (23) 

Following Ref. [38l . a higher-derivative sum rule can be 
obtained by multiplying S~ by iS^~ and expressing the 
expectation value by Eq. (|22[) , 

q i=l,2 

= 2 £ Jn (c» - i<S 2 )) . (24) 

Considering the ground state, at T = 0, we have the 
exact results 

C-+(0) = 0, C£*(0) = i (^)(0) = i (25) 

Here C~ + (0) is independent of n, so A g i j2 (0) must be 
independent of q. This requires K„ tm (0) = which 
yields 0^2(0) = 1- Concerning the zero-temperature 
values of a\ z 2 , they can be determined only in the limit 
T -» since' Eqs. © and © for C^ z contain with 



limr-i-o A^g Z = 0. To evaluate the thermodynamic prop- 
erties, the correlators C* z and C~ + for n ^ are calcu- 
lated by Eqs. © and ([22]). whereas (5 12 ) and the vertex 
parameters a^J an( i a i Z 2 are determined from the sum 
rules [Eqs. dMJ), and = 1/4], the equality © 

and the regularity condition (|21l) . 



III. RESULTS AND DISCUSSION 

In what follows we put | Ji | = 1 . The coupled system of 
nonlinear algebraic self-consistency equations derived in 
Sec. HJis solved numerically using Broydcn's method^ 
which yields the solutions with a relative error of about 
10~ 7 on the average. The momentum sums for N — > 
00 are transformed to integrals which are evaluated by 
Gaussian integration. Tracing the GFM solution to very 
low temperature we find it to become less trustworthy 
for J2 approaching J2 = 0.25. Therefore, below we will 
present GFM results for J2 < 0.2 only. On the other 
hand, there is no restriction with respect to the value of 
J r 2 for the ED calculations. 

We restrict our discussion to magnetic fields h < 0.1. 
For typical values of the exchange parameters of the order 
of 10 me V 20 ' 21 that corresponds to values of the magnetic 
field H < 10T mostly used in experiments. 



A. Magnetization 

First we consider the magnetization (S z ), see Fig. [1] 
Obviously, at fixed magnetic field the frustration reduces 
the magnetization for T > 0, whereas an increase of the 
magnetic field at fixed J2 leads to an increase of (S z ). 
The GFM and ED data agree well with each other. 



B. Magnetic susceptibility 

The uniform magnetic susceptibility \ = d(S z )/dh at 
zero field diverges at T — indicating the ferromagnetic 
phase transition. Based on GFM data it was found in 
Ref. 03 that lim X T 2 = (l-4J 2 )/24. (Note that for J 2 = 

the rigorous Bethe-ansatz result^ is reproduced.) In 
finite magnetic fields one has lim v = 0, and y exhibits 

T->0 

a maximum. The position T^ ax and the height Xmax of 
this maximum depend on the strength of the magnetic 
field h and also on the frustration J2. 

The field dependence of the position of the susceptibil- 
ity maximum was discussed, for example, in connection 
with experiments on Lao. giMno. 95031^ For the unfrus- 
trated ferromagnet (J2 = 0) it was show n 37 ' 38 ' 46 that 
interestingly power laws for the field dependence of T^ ax 
and Xmax are valid, namely 

T max = ah? , Xmax = bh e . (26) 
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FIG. 2: Magnetic susceptibility \ versus temperature T 
(GFM - lines; ED N = 20 - filled circles; ED N = 16 - 
open triangles) for (a) J2 = 0.05 and (b) J2 = 0.15 with 
h — 0.01, 0.05, 0.08, from top to bottom. The insets show the 
position and height of the maximum obtained by GFM (filled 
squares) and ED N — 20 (open circles) as well as the fit by 
the power law (lines), see Eq. (|26l) . 



FIG. 3: Field dependence of (a) the position T£ ax and (b) the 
height x-max of the maximum of the magnetic susceptibility 
obtained by GFM (main panels). The lines in the upper panel 
show the fits by the power law Eq. (|26|) . In the insets the 
GFM results (filled triangles) and the corresponding ED data 
for N = 20 (open triangles) are compared for J2 = 0.15. 



We mention that our analysis of the susceptibility data 
reproduces the results reported in Refs.l37l38 for J 2 = 0. 
Note that a, 7, b, and /3 depend on the spin quantum 
number s and dimensionality Di 37 i 38 

Now we discuss the susceptibility for the frustrated 
ferromagnet, see Fig. [5J The behavior is similar to the 
unfrustrated case, i.e. for a fixed J 2 , the position T* ax 
of the maximum increases and the height Xmax decreases 
with increasing field. For a fixed h with increasing of 
frustration J 2 , the position of the maximum is slightly 
shifted to lower temperatures, whereas the height is al- 
most fixed for J 2 < 0.2. The ED and GFM data match 
to each other reasonably well except the noticeable dif- 
ference near the maximum. 

Next we discuss the validity of the power law (j2"6")) for 
finite frustration J 2 > 0, see insets of Fig. [2] We find that 
the GFM as well as the ED data are well described by 
the power law ([26]) in the whole range of frustration (i.e. 
J 2 < 0.2 for the GFM and J 2 < 0.25 for the ED), see the 
log-log plots of T* ax and Xmax in Fig. [3J The influence 
of frustration is quite weak. To obtain the dependences 



of the coefficients on the frustration J 2 we fit the position 
and the height of the maximum by the power law using 
the data points shown in Fig. [3] The variation of the 
coefficients a, 7, 6, and /3 with J 2 shown in Figs. [4] and 
[5] is weak for J 2 < 0.2 (note the enlarged scale in Figs.0] 
and [5]). The difference between the ED and GFM data 
may be attributed to finite-size effects. Obviously, the 
ED for N = 20 match better to the GFM data than 
the ED data for N = 16 (except for the parameter a at 
larger J 2 .) For comparison, we have also shown the exact 
Bethe-ansatz results 3 -! for J 2 = 0. The GFM data agree 
well with the exact data, except for the parameter a, 
where the difference is about 23%. For J 2 > 0.2, we have 
only ED data. In this parameter region the finite size- 
effects are very small, and there is a noticeable change of 
the coefficients a, 7, and j3 with frustration. 



C. Correlation functions 

Next we consider the spin-spin correlation functions 
depicted in Fig. [5] Obviously, the GFM results for the 
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FIG. 4: Coefficients (a) 7 and (b) a of the power law T^ ax = 
ah 1 in dependence on the frustration J2 obtained by GFM 
and ED (N=16,20). For J 2 = also the Bethe-ansatz result 
is shown. 



FIG. 5: Coefficients (a) /3 and (b) 6 of the power law Xmax = 
bh 13 in dependence on the frustration J2 obtained by GFM and 
ED (N=16,20). For J 2 = also the Bethe-ansatz result is 
shown. 



correlation functions agree quite well with the ED data. 
Due to the magnetic field in z-direction we have a dif- 
ferent behavior of the longitudinal correlation functions 
(S^S 2 ) and the transverse ones (SgS%) = (S#S»>. At 
T ^> h this difference disappears. At low temperatures 
the transverse correlation functions exhibit a maximum, 
where its position and height depend on the spin-spin 
separation n, the frustration J 2 , and the magnetic field 

h. As expected, frustration suppresses the correlations, 

i. e., there is a faster decay of the correlation functions 
with increasing J 2 . The strength of the magnetic field 
at fixed J 2 influences the correlation functions at low 
temperatures only. Increasing of h yields a shift of the 
maximum in (SqS%) to higher temperatures and a slower 
decay of the correlation functions with T. 



D. Correlation length 

Within the GFM formalism we are able to calculate 
the correlation length i 10 i 38 ' 39 Due to the field-induced 
anisotropy we have again to distinguish between the lon- 
gitudinal correlation length ^ and the transverse cor- 
relation length £ 22 . Both quantities can be obtained by 



expanding the static susceptibilities 3 ^ Xq^t v ^ = "I — i zz i 
around the magnetic wave vector q = which leads to 



x?" = xo7{i + (€"") V}. 



We obtain 



■Jid- 4J 2 C 2 A" 1 



— 



(S')h 

= E 

n,m,(=l,2) 



2h 2 



n A 1 J K 



+ - 



and 



2(Ji+4J 2 )(j 1 afCf 



4J 2 af 



Co 



(27) 
(28) 

(29) 



where A 22 is given by Eq. ©. As found in Ref. [H, the 
transverse and longitudinal correlation lengths have qual- 
itatively different temperature dependences at J 2 = 0. In 
particular, the longitudinal correlation length £ 22 reveals 
an anomaly, namely a shoulder at T ~ 0.4 appearing at 
low magnetic fields, cf. also Fig. [Tfb). 

Considering the transverse correlation length £ H 
shown in Fig. [TJa), the magnetic field cuts off the di- 
vergence of the zero-field correlation length at T = 
(cf. Ref. 01). This is related to the absence of the 
T = phase transition if h > 0. From Eq. (|27l) we 
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0.25 




FIG. 6: Spin-spin correlation functions (So^n) and (SES%) 
for n = 1 and n = 4 for (a) h = 0.03 and (b) ft = 0.1 (GFM 
- lines, ED data for N — 20 - symbols, black lines/symbols: 
Ja = 0, n = 1; red lines/symbols: J2 = 0.15, n = 1; blue 
lines /symbols: J2 = 0, n = 4; green lines/symbols: J2 = 0.15, 
n = 4). 



FIG. 7: (a) Transverse and (b) longitudinal correlation 
lengths in dependence on the temperature for frustration pa- 
rameters J2 = 0, 0.1, 0.15 (from top to bottom) and for various 
values of the magnetic field. 



E. Specific heat 



find £ + ~(T = 0) = ^/ °- 5 ~ 2j2 , i.e., at low temperatures 

the transverse correlation length increases as h~ x / 2 for 
h — > 0, and it decreases with increasing frustration ac- 
cording to (1 — 4J2) 1 / 2 . At temperatures larger than 
the magnetic field the dependence on the field strength 
becomes weak. As already discussed for the correlation 
functions, at those temperatures both correlation lengths 
£ H and £ zz approach each other indicating that the 
magnetic field becomes irrelevant. On the other hand, 
the decrease in and £ 22 with increasing frustration 
is observed in the whole temperature range shown. At 
low temperatures T < h we find £ 22 < ^ . At very 
low temperatures this relation might be reversed. Un- 
fortunately, the GFM does not allow here a conclusive 
statement on £ 22 as T — > 0, since in Eq. (|29| for £ 22 both 
the nominator and denominator go to zero for T — > 
leading to numerical uncertainties. 



Let us first recapitulate some relevant results on the 
specific heat C(T) = du/dT with u = (H)/N = 
Ji (Cf + + Cf z ) + J 2 (C7+ + Cf 2 ) -h(S z ) from previous 
investigations ^il2^L^ For the unfrustrated ID s = 1/2 
ferromagnet at h = the specific heat shows a broad 
maximum at T ~ 0.35, see, e.g., Refs. [i"oH37l and[f| In a 
weak magnetic field h < h* a double-maximum structure 
appears. For the field h* below which the additional low- 
temperature maximum appears, a Bethe-ansatz analysis 
gives h* = 0.008, see Ref. [13, whereas the GFM yields 
h* fn 0.07>2L22 On the other hand, for the J x -J 2 model 
in zero field also a double-maximum structure in C(T) 
was found. In this case it is induced by frustration, and 
it appears within the GFM for J 2 > 0.16f^ Moreover, it 
was found in Ref. [l^ that, although the existence of an 
additional low-temperature maximum can be detected by 
GFM and ED, its position and height are strongly size 
dependent, if the maximum is located at very low tem- 
peratures. 

Our results for the specific heat of the frustrated fer- 
romagnet in a magnetic field are shown in Fig. [5] Analo- 
gously to Refs . l37l and l38l. a double- maximum structure is 
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FIG. 8: Specific heat for (a) J 2 = 0.05 and (b) J 2 = 0.15 
with h = 0.005, 0.03, 0.05, 0.08, 0.1, from bottom to top (lines 
- GFM; symbols - ED for N = 20). 



found for small field strengths. An increase of the mag- 
netic field shifts the additional low-temperature maxi- 
mum to higher temperatures and increases its height. 
At a sufficiently high field strength ft > ft* the typical 
broad maximum of the pure ferromagnet at zero field 
vanishes, since the low-temperature maximum spreads 
into the region of the former broad maximum. Although 
the GFM and the ED data for the specific heat differ 
quantitatively, the above mentioned qualitative behav- 
ior is found by both methods., cf. also the discussion in 
Refs. Holl37l and HH Recalling that for the unfrustrated 
ferromagnet the double- maximum structure disappears 
for ft > ft* « 0.07 (ft > ft* w 0.008) calculated by GFM 
(ED), cf. Ref. [H it can be seen in Fig. Mb) that ft* is 
shifted to higher values due to frustration. Thus, in Fig. 
MJ°) at J2 = 0.15, the double-maximum structure is still 
visible in the GFM data at h = 0.08 and in the ED data 
at ft = 0.03. 

As pointed out in Refs. and H^, for J2 = the 
position and the height of the additional low-temperature 
maximum follow a power law 



T- 



.1 = rh" , 



dh s . 



(30) 



To determine the coefficients r, p, d, and S in dependence 
on the frustration J2 we follow Ref. [H and fit the po- 



FIG. 9: (a) Coefficients p and r of the power law T^ ax ^ = rh p 
in dependence on the frustration parameter J2 obtained by 
GFM. For comparison the Bethe-ansatz results are shown for 
J2 = 0. (b) Coefficients S and d of the power law C ma x,i = 
dh s in dependence on the frustration parameter J2 obtained 
by GFM. For comparison the Bethe-ansatz results are shown 
for J 2 = 0. 



sition and the height of the low-temperature maximum 
at low fields, h = 0.001 - 0.01 in steps of 0.001, by the 
power laws of Eq. ([3H]) . For these values of the magnetic 
field, the ED maximum in C(T) quite strongly depends 
on the size N. Therefore, in Fig. Owe present the GFM 
results only. It is obvious that up to J 2 ~ 0.15 all coeffi- 
cients depend weakly on J2. Only beyond J2 ~ 0.15 the 
coefficients r and p increase noticeably. 

Let us finally comment the obvious differences in the 
GFM and ED data for the specific heat, see Fig. [SI The 



10.36-41 



is a 



GFM used here and in many previous papers 
quite universal analytical method which, as a rule, can- 
not yield quantitatively highly accurate data. On the 
other hand, the accuracy of the ED data is also limited 
at low temperatures due to significant finite-size effects, 
see, e.g., the discussion in Ref. [l^- In cases where other 
accurate methods (for instance the quantum Monte Carlo 
method which is applicable to unfrustrated quantum spin 
systems, see e.g. Ref. l38l ) are available, it was found 
that GFM results for the specific heat at low temper- 
atures show quite large quantitative differences to pre- 
cise data, whereas other quantities, such as the suscep- 
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tibility, are in better agreement. However, it has been 
demonstrated that both the GFM and the ED in gen- 
eral yield qualitatively correct results also for the spe- 
cific heati 10 ' 38 In particular, the existence of the addi- 
tional low-temperature maximum in the specific heat at 
weak magnetic fields discussed in this section was con- 
firmed by a Bethc-ansatz analysis 3 - 7 - as well as a quan- 
tum Monte Carlo calculation 3 ^ for the unfrustrated ID 
s=l/2 ferromagnct, where both methods are applicable. 
Hence, we argue that the existence of the additional low- 
temperature maximum in the specific heat is not ques- 
tioned, however, quanitative statements have to be taken 
with caution. 



IV. SUMMARY 

In this paper we use exact diagonalization of finite 
chains and a second-order Green's function method to 
study the influence of a frustrating next-nearest-neighbor 
coupling Ji on low-temperature thermodynamic quanti- 
ties, such as magnetization, correlation functions, trans- 
verse and longitudinal correlation lengths, spin suscepti- 



bility, and specific heat of the ID s = 1/2 J1-J2 Hcisen- 
berg ferromagnct in a magnetic field h. We consider 
J2 < 0.25|Ji|, i.e., the ground state of the model is fer- 
romagnetic, but the low-temperature thermodynamics is 
strongly influenced by the frustrating J2. The results of 
both methods are in good overall agreement. 

The position T* ax and the height Xmax of the max- 
imum of the uniform susceptibility appearing in finite 
magnetic fields can be described by the power laws 
Tmax = a ^ 7 an d Xmax = bh@ , where the coefficients a, 
7, b, and (3 weakly depend on Ji. The double-maximum 
structure of the specific heat C(T) found for J2 = at 
low field a 37 ' 38 occurs also for J2 7^ 0. The field strength 
at which this structure vanishes is slightly increasing 
with frustration. The position and the height of the 
additional low-field low-temperature maximum in C(T) 
in dependence on the field strength can be also described 
by power laws. 
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